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ABSTRACT 

We have made a likelihood statistical analysis of the angular correlations in the COBE- 
DMR two-year sky maps by Monte Carlo simulation of the temperature fluctuations. 
We assume an open universe and consider as primordial power spectrum the Harrison- 
Zeldovich one, P{k) = Ak. We hud that the flatness of the universe is not implied by 
the data. The quadrupole normalization amplitude, Qrms-PS, is related to the density 
parameter, SI, by Qrms-PS = 10.67 -|- 55.81S1 — 128.59S1^ + 81.26S1^ fiK. We have deter¬ 
mined the p.d.f. of SI due to cosmic plus sampling (i.e. 20° galactic cut) variance which 
generically shows a bimodal shape. The uncertainty as given by the r.m.s. is ss 0.35, 
therefore to better constrain SI experiments sensitive to higher multipoles (/ > 20) 
should be considered. 
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1 INTRODUCTION 

Rapid growth of the number of cosmological observations re¬ 
veals the week points of standard Q = f CDM models. These 
models are unable to simultaneously account for the follow¬ 
ing observations: galaxy and cluster autocorrelation func¬ 
tions, large scale streaming motions and the normalization 
amplitude of the Cosmic Microwave Background (CMB) 
spectrum. A low-density (Q < f) universe is favoured by 
many of the observations (see, e.g.. Coles & Ellis 1994). 
Other possibilities are the flat universe with a cosmologi¬ 
cal constant, tilted models and mixed hot-cold dark matter 
models. Flat and low density models present quite differ¬ 
ent shapes of the CMB anisotropies due to geometric effects 
(Sugiyama & Silk 1994). We can use the angular correlation 
of CMB maps to test open, flat and low density flat models. 
By assuming a Harrison-Zeldovich (HZ) spectrum, Bunn & 
Sugiyama (1995) constrained the value of the cosmological 
constant from COBE-DMA maps. 

Open universes have not been studied in detail because 
the precise shape of the density power spectrum beyond the 
curvature scale in these universes has not been determined 
dehnitively. The pioneering calculation of Wilson (1983) as¬ 
sumed a power law in wavenumber k which is dehned as 
— 1 where q is the eigenvalue of the Laplacian, while 
Kamionkowski and Spergel (1994) considered power laws in 
spatial volume and in scalar spatial Laplacian eigenvalue. 
Recently, the power spectrum has been calculated by as¬ 
suming inflation at the early epoch of the universe (Lyth & 


Stwart 1990; Ratra & Peebles 1994). This model has been 
confronted with some observations which do not include all 
of the COBE-DMB experimental features (Kamionkowski et 
al. 1994; Sugiyama and Silk 1994) but Gorski et al. (1995), 
incorporating such features, conclude that the COBE-DMB 
two-year data do not provide sufficient statistical power to 
prefer a particular value of Q. However, further modihcation 
of the power spectrum might be necessary because of the 
bubble nucleation in the early epoch which allows the exis¬ 
tence of the open universe together with inflation (Bucher, 
Goldhaber & Turrok 1995; Lyth & Woszczyna 1995; Ya¬ 
mamoto, Sasaki & Tanaka 1995). 

In this paper, we study the angular correlation of CMB 
anisotropies with the simplest possible shape of the initial 
power spectrum, i.e., the HZ spectrum P{k) ex. k where 
k is the wave number. CMB large scale anisotropies in a 
flat universe are produced by density fluctuations on the 
last scattering surface (Sachs & Wolfe 1967). At angular 
scales greater than {2Q}P)° a generalized formula for curved 
spaces can be obtained (Anile and Motta 1976). In this case 
the anisotropies can be expressed in terms of gravitational 
potential fluctuations at recombination and an integrated 
effect of the time-varying gravitational potential along the 
photon trajectory that depends on curvature (Abbott and 
Schaefer 1986, Traschen and Eardfey 1986, Gouda et af. 
1991). In this work we have calculated CMB anisotropies by 
numerically integrating the generalized formulas for curved 
spaces (Sugiyama & Silk 1994) as shown in Fig. f. Tak¬ 
ing a muftipofe expansion for the temperature fluctuations 
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Figure 1. The angular power spectrum for various Q values: 1 
(solid), 0.8, 0.6 (upper and lower dashed lines), 0.4 (dash-dotted), 
0.3, 0.2, 0.1 (lower, middle and upper dotted lines. The Harrison- 
Zeldovich primordial spectrum was assumed. 


AT/T = E l,m 0^1,171 Yim, only niultipoles up to 30 are sensi¬ 
tive to the COBE-DMB. angular resolution. 

Analysis of the hrst year COBE-DMB. maps using an¬ 
gular correlations resulted in a spectral index, ra of 1.1 ±0.5 
(assuming a power law spectrum P{k) = Ak"), and a rms- 
quadrupole-normalized amplitude, Qrms-ps, of 16 ± AfiK 
(Smoot et al.l992). An analysis based on the ‘genus’, or to¬ 
tal contour curvature of anisotropy spots, found n = 1.2±0.3 
(Torres 1994a, 1994b). This result was obtained by htting 
the coherence angle of the COBE-DMB temperature maps 
for a hxed Qrms-PS of 16fiK, assuming a power law primor¬ 
dial density fluctuation spectrum. By allowing variations 
of the Qrms-PS (oc parameter, a relation between 

Qrms-PS and n was also found (Smoot et al. 1994, Torres et 
al. 1995). Analysis of the hrst two years COBE-DMB maps 
shows that the Harrison-Zeldovich spectrum remains consis¬ 
tent with the spectral index predicted by the data (Bennett 
et al. 1994, Gorski et al. 1994). In all the previous analyses 
a flat universe with Q = 1 was assumed. 

Even in an ideal noise-less experiment the measured 
parameters will have an uncertainty from cosmic variance. 
We have used Monte Carlo simulations to study the extent 
to which cosmic and sampling variance obscures the infor¬ 
mation about Q which can be extracted from a correlation 
analysis. 

In section 2 we present the COBE-DMR Monte Carlo 
simulation methodology while in section 3, we obtain the 
minimum range of Q implied by cosmic and sampling vari¬ 
ance. In section 4 we apply the correlation technique to the 
COBE-DMB data for low-Q models and in section 5 sum¬ 
marize and discuss the results. 


2 MONTE CARLO METHODOLOGY 

Taking into account instrument noise, non-uniform sky cov¬ 
erage, galactic cut, smearing, pixellization scheme, and the 
DMR beam characteristics we simulated a set of cosmologi¬ 
cal models dehned by Q, Qrms-PS, and HZ primordial spec¬ 
trum. We generated 3200 simulated CMB sky maps using a 
harmonic expansion of the temperature for each of the 6,444 


DMR pixels: 

i 

/\T(e,<j>) = EE km\bi,m cos(m(j)) 

1 = 2 m = 0 

+bi-m sm{m(l>)]NrWiPricos9), (f) 

where 

r (2l±f)(l-m)! l^^^^ 

* 47r(l ± m)! 

km — t/ 2 for m 0 and ko = 1. PE {cos 9) are the associated 
Legendre polynomials. The coefficients bi^rn are real stochas¬ 
tic, Gaussianly distributed variables with zero mean and 
model dependent variance {bf .^} that are easily obtained 
from the multipole coefficients Ci {{bi .^} = CiEir) given 
by Sugiyama & Sifk (4994) in their Figure fa. The weights, 
W{, for DMR given by Wright et af. (4994) are used and the 
series is truncated at I = 30. A smaff beam smearing cor¬ 
rection is appfied which accounts for the spacecraft motion 
during the i/2 second integration time. For each reafization 
two maps are generated by adding to the cosmic signaf the 
noise corresponding to the channef combination (A±B)/2 of 
COBE-DMR 53 GHz and 90 GHz. The noise is determined 
by instrument sensitivity and the number of observations 
per pixef. The angufar cross-correfation of the 53 GHz and 
90 GHz maps for |6| > 20'’ is cafcufated and binned in 36, 
5'’ bins in the manner of the COBE-DMR data (Bennett et 
af. 4994). A combination of the two independet frequency 
maps hefps reducing the introduction of systematic effects. 

The standard fikefihood statistic L associated to the 
cross-corefation is (ra = 36) 


E = 


(27r)"/2(detM)i/2 


(2) 


36 36 

= EE(<^') - (3) 

t=i j=i 


(Ct) is the average of the cross-correfation for the 3200 re- 
afizations at bin i. Cp is the cross-correfation for the 
COBE-DMR maps at angufar scafe i (note that G must 
not be confused with the muftipofe coefficients Ci). Mtj is 
the covariance matrix cafcufated with the Monte Garfo re- 
afizations: 


- Nr 


(Cf-(G))(C''-(C,)). 


( 4 ) 


The fikefihood method is preferabfe to the one 
based on the foifowing test: we have used as ‘data’ 400 
reafizations from each of the foifowing Q — Qrms-ps{p^) 
pair of vafues: O.f - 45,0.2 - 47,0.3 - 49,0.4 - 47,0.6 - 
45,0.8 — i5,i — i9 (these are the modefs with maximum L 
for hxed Q, see section 4). For each Q we get the distri¬ 
bution of Q-vafues for which L and nre maximum and 
minimum, respectivefy, by comparing each one of the 400 
‘data’ reafizations with 3200 reafizations of the modefs with 
Qrms-PS = 7, i2,i5,i7, 49,24,23,25,27. The mean vafue 
of the L method is within 0.5 units of the true vafue with a 
dispersion of aq Si i.5 for aff Q vafues. However the mean 
vafue of the x^ method systematicaffy overestimates the true 
Qrms-PS vafue in Si 3 units and with a signihcantfy bigger 
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dispersion. So, the likelihood method accurately predicts the 
Qrms-ps value with a relatively small uncertainty whereas 
the has systematic errors. This result could be explained 
by the fact that the number of bins (i.e. 36) begins to be 
high enough for the likelihood to be unbiased. 


3 UNCERTAINTY FROM COSMIC PLUS 
SAMPLING VARIANCE 

Assuming a noise-less radiometer, we now study the cosmic 
variance associated to the cosmological parameters derived 
from the correlation analysis. The COBE-DMB. radiometer 
noise is taken into account in the next section when applying 
the method to the real data. 

The non-ergodic character of the CMB temperature 
random held on the cosmic sphere limits the accuracy in 
determining the statistical properties of the held from a sin¬ 
gle realization (i.e. averages over the celestial sphere do not 
coincide with statistical averages over the ensemble of real¬ 
izations). Even in the ideal case when the cosmic background 
temperature is measured in every part of the celestial sphere 
with no galactic cut and with negligible noise we would only 
be able to measure the expected values of the held within 
certain errors, the cosmic variance. Uncertainties in the two- 
point correlation function and in other higher moments due 
to cosmic variance have been calculated (Scaramella and 
Vittorio 1990, Cayon, Martinez-Gonzalez and Sanz 1991, 
White, Krauss and Silk 1993, Srednicki 1993, Gutierrez de 
la Gruz et al. 1994). Uncertainty in the determination of the 
spectral index n for a flat universe has been also estimated 
in the genus analysis (Torres et al. 1995). 

We calculate the cosmic plus sampling variance uncer¬ 
tainty when analysing the cross-correlation of two hypothet¬ 
ical sky maps with the characteristics of the COBE-DMB. 
experiment (i.e. 7° FWHM beam) assuming a noise-less ra¬ 
diometer and |6| > 20° galactic cut. Here we have two free 
parameters, Q and Qrms-PS, and seven Monte Garlo data 
sets were generated with 3200 realizations per set, one for 
each Q value - 0.1, 0.2, 0.3, 0.4, 0.6, 0.8,and 1.0 - , with hxed 
Qrms-PS since the results are not too sensitive to changes 
in this parameter over a realistic range of values. To esti¬ 
mate the uncertainty in Q due to cosmic variance we com¬ 
puted the likelihood L for each of 400 realizations from a 
previously selected model. As a hrst case Q and Qrms-PS 
were taken to be 0.4 and 17 fiK (this is the normalization 
for which L is maximum for the model Q = 0.4, as ob¬ 
tained in the next section). Cf * in the likelihood dehni- 
tion was replaced by Cf, the value of the cross-correlation of 
the A;-th realization in bin i. For each realization the value 
of Q which maximized L (fim) was found. The resulting 
distribution of Qm is: 0, 11, 11, 44, 0, 0, 34% probabilities for 
Q = 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1. Therefore, the flat model is 
indistinguishable from the ‘true’ 0.4 model and only models 
with Q = 0.1, 0.6, 0.8 are clearly distinguished. This counter¬ 
intuitive result can be accounted for by the similarities in the 
angular power spectrum shown by the Q = 0.4 and 1 models 
and the large differencies between Q = 0.4 and 0.8 models 
(see Fig. 1). In general, for other values of the ‘true’ Q a 
similar bimodal distribution appears. 


Figure 2. The likelihood function L (normalized to 1 at the max¬ 
imum value) as a function of Q and Qrms — PS derived from the 
study of the cross-correlation of the of the 53 GHz and 90 GHz 
maps for |6| > 20°. We have performed a bilinear interpolation 
in the {Qrms — PS: U) plane. The contours are 10 sections at 0.1 
intervals of the maximum L-value. 


4 ANALYSIS OF THE GOBF-DMR DATA 

We used two-year data from the most sensitive radiometers 
at 53 and 90 GHz. To determine the limits on Qrms-PS and 
Q imposed by the COBE-DMB cross-correlation, a grid of 
Monte Garlo data sets was generated for seven selected Q 
values (0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1) and Qrms-PS between 
7 and 30 /tK in 2 /tK steps. For each simulated (A-|-B)/2 sky 
map at 53 and 90 GHz, the cross-correlation was computed 
and the L statistic calculated. The number of realizations 
was 3200, a sufficient number as found after testing for the 
convergence of the relevant quantities. 

Given the Harrison-Zeldovich primordial spectrum, we 
searched for the model (Qrms-PS, fl) that maximized the 
likelihood L using Monte Garlo data sets. We hud a degen¬ 
eracy in the plane (Qrms-PS,^) that maximizes the value 
of L. To obtain such a relation between Qrms-PS and Q, 
an error to the value of Qrms-PS with maximum L for 
hxed Q is assigned according to the spread of Qrms-PS val¬ 
ues obtained from the likelihood method with each realiza¬ 
tion from that Qrms-PS model taken as the input data. 
We then ht a cubic polynomial to the pairs (Qrms-PS,^) 
with maximum L considering the corresponding error bars. 
A qualitative origin of such a relation can be understood 
in terms of the combined contributions of potential fluc¬ 
tuations and the integrated effect due to curvature which 
results in a partial cancellation near Q Si 0.4. The re¬ 
lationship between our estimates of Q and Qrms-PS can 
be approximated by a cubic polynomial: Qrms-PS (Q) = 
10.67 -|- 55.81Q — 128.59Q^ -|- 81.26Q^ /tK. The maximum 
L is obtained for Q = 0.1, Qrms-PS = 15/tK, although the 
likelihood is not signihcantly worse at higher values of Q 
or near 1. In fact, all {tl,Qrms-p,s) values satisfying the 
previous cubic relationship are statistically similar. Figure 
2 shows L versus Q and Qrms-PS and also 10 contours at 
levels 0.1 of the maximum L-value. Assuming a Bayesian 
interpretation of L, one can estimate the 2D regions of con- 
hdence at different levels (la, 2a, and Ocr) as seen in Fig. 3 
where it is clear that above 2a all Q models are statistically 
equivalent. Thus, GOBF-DMR two-year data do not favour 
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Figure 3. Confidence regions at Ic, 2(T and 3(T obtained from the 
likelihood L assuming a Bayesian interpretation. Above the 2(T 
level, there is no any favoured value of Q. 

a flat universe with Q = 1. Taking into account an analo¬ 
gous result by Gorski et al (1995) for a different primordial 
power spectrum (the Ratra-Peebles one), in the sense that 
the COBE-DMB. two-year data are not able to provide any 
specihc answer about Q, we can also conclude that the men¬ 
tioned result is not specihc of the HZ power spectrum but 
can be considered as a general result independently of the 
primordial power spectrum assumed. 

In table 1 we give the relationship between the nor¬ 
malization amplitude A of the primordial power spectrum 
and the values of the parameter Qrms-ps obtained in our 
maximum likelihood analysis. The corresponding rms mass 
huctuation at 8h~^ Mpc, as, for the different Q models con¬ 
sidered here is derived for several pairs of values of Qb, h. It 
is clear that all models with Q < 0.3 are inconsistent with 
observations of clustering at small scales. 

Table 1. Normalization of the primordial power spectrum A for 
different values of Q. A(20/iK) is the amplitude in P{k) = Ak 
that gives an rms quadrupole amplitude of 20/iK, while A[DMR) 
refers to the COBE-DMR quadrupole normalization. A is given in 
units Mpc^. The resulting eg (h, Qb) is also given. Qrms — PS 
is the value favoured by our analysis for the different values of Q . 


5 CONCLUSIONS 

Based on the standard likelihood statistic, we have consid¬ 
ered the test for the temperature cross-correlation of sky 
maps and we are able to constrain the quadrupole normal¬ 
ization amplitude Qrms-PS and the density parameter Q 
of the density fluctuation power spectrum at recombination 
with the COBE-DMR two-year data. It is found that the 
two parameters lie within a region around the polynomial 
Qrms — PS (Q) = 10.67-|-55.81Q- 128.59Q^ -|-81.26Q^ fiK for 
a Harrison-Zeldovich spectrum . 

We have also studied the effect of cosmic plus sampling 
variance on the determination of the density parameter Q 
when using the cross-correlation as the statistical quantity 


in the comparison with the COBE-DMR data. The result 
depends on the ‘true’ Q model assumed but generically a 
bimodal distribution appears, i.e. there exists another Q- 
model which is indistinguishable from the ‘true’ one. 

We can also conclude that the scale invariant spectrum 
predicted by inflation and the flatness of the universe are 
neither ruled out nor directly implied by the COBE-DMR 
two-year data. 
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Slo 

T(20//K) 

(xlO®) 

A(DMR) 

(xlO®) 

(78(0.5,0.03) 

(78(0.5,0.06) 

(78(0.8,0.03) 

(78(0.8,0.06) 

Qrms — PS 

ifiK) 

0.1 

2.894 

1.638 

0.0740 

0.04245 

0.1283 

0.06493 

15.05 

0.2 

1.590 

1.195 

0.2037 

0.1578 

0.3580 

0.2641 

17.34 

0.3 

1.334 

1.084 

0.3566 

0.3011 

0.6185 

0.5061 

18.03 

0.4 

1.353 

1.050 

0.5250 

0.4620 

0.8952 

0.7703 

17.62 

0.6 

1.894 

1.125 

0.9197 

0.8433 

0.1520 

1.375 

15.42 

0.8 

2.179 

1.165 

1.318 

1.233 

0.2121 

1.965 

14.63 

1.0 

0.798 

0.732 

1.337 

1.268 

0.2086 

1.965 

19.15 
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